Spontaneous ferromagnetism in the spinor Bose gas with Rashba spin-orbit coupling 
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We show that in the two- component Bose gas with Rashba spin-orbit coupUng an arbitrarily small 
attractive interaction between bosons with opposite spin induces spontaneous ferromagnetism below 
a finite critical temperature Tc. In the ferromagnetic phase the single-particle spectrum exhibits a 
unique minimum in momentum space in the direction of the magnetization. For sufficiently small 
temperatures below Tc the bosons eventually condense into the unique state at the bottom of the 
spectrum, forming a ferromagnetic Bose-Einstein condensate. 
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I. INTRODUCTION 

Due to recent progress in the field of ultracold gases, 
spinor Bose-Einstein condensates with various types of 
spin-orbit coupling can now be realized experimentall}^^ ^ 
by using spatially varying laser fields to couple in- 
ternal pseudo-spin degrees of freedom to the momen- 
tum. These experiments have motivated many theoreti- 
cal investigations of spin-orbit coupled multi-component 
Bose gases. ^"^^ Of particular interest have been two- 
component bosons with isotropic Rashba-type^^ spin- 
orbit coupling, where in the absence of interactions the 
energy dispersion assumes a minimum on a circle in mo- 
mentum space. ^ If the bosons do not condense, such a 
surface in momentum space can be called a Bose sur- 
face, in analogy with the Fermi surface of an electronic 
system.E^EH Of course, bosons do not obey the Pauli ex- 
clusion principle, so that the Bose surface cannot be the 
boundary between occupied and unoccupied states; how- 
ever, the Bose surface defines the location of the low- 
energy excitations in the system, similar to the Fermi 
surface of an electronic system. 

The fact that Bose-Einstein condensation (BEC) in 
Bose systems where the dispersion has degenerate min- 
ima on a surface in momentum space differs qualitatively 
from conventional BEC has been pointed out a long time 
ago by Yukalov.^^ He studied BEC in an interacting Bose 
system whose energy dispersion is minimal on a sphere in 
momentum space. Assuming that the bosons condense 
with equal weight into all states on this sphere, he showed 
that the condensed state does not exhibit off-diagonal 
long-range order and is also not superfluid. 

Due to the degeneracy of the single-particle energy in 
the spinor Bose gas with Rashba-type spin-orbit cou- 
pling, in the non-interacting limit BEC is prohibited at 
any finite temperature. Hence, finite temperature BEC 
in this system must be an interaction effect .1^^ Interac- 
tions are expected to remove the ground state degener- 
acy and various types of exotic ground states have been 
proposed. "^J^ Which phase is realized experimentally de- 
pends on the specific properties of the interaction. At 



this point a generally accepted agreement on the nature 
of the ground state in the spinor Bose gas with Rashba- 
type spin-orbit coupling has not been reached. 



Phase transitions in systems whose fluctuation spec- 
trum exhibits minima on a surface in momentum space 
form their own universality class, the so-called Brazovskii 
universality class p^ for example, the critical behavior in 
cholesteric liquid crystals belongs to this class .122 Because 
scaling transformations and mode elimination in renor- 
malization group calculations should be defined relative 
to the low-energy manifold, the classification of interac- 
tion vertices in systems belonging to the Brazovskii uni- 
versality class is different from the corresponding classi- 
fication in systems where the low-energy manifold con- 
sists of a single point; in particular, all two-body scat- 
tering processes where the momenta of the incoming and 
the outgoing particles lie on the low-energy manifold are 
marginal, so that in renormalization group calculations 
one should keep track of infinitely many marginal cou- 
plings. Note that also normal fermions belong to the 
Brazovskii universality class, because the Fermi surface 
can be identified with the low-energy manifold relative 
to which scaling transformations should be defined.^"^ 
Two-component bosons with Rashba-type spin-orbit cou- 
pling are therefore another example for a quantum sys- 
tem which belongs to the Brazovskii universality class. 



In this work we shall further investigate interaction ef- 
fects on spinor Bose gases with Rashba-type spin-orbit 
coupling. We shall consider the specific case where the 
interaction g^ between two bosons with opposite pseudo- 
spin is attractive. We find that in this case for any finite 
density an arbitrarily small interaction ^± < leads to 
spontaneous ferromagnetism below some finite tempera- 
ture Tc. In the ferromagnetic phase, the single-particle 
dispersion has a unique minimum in momentum space, 
so that the bosons eventually condense at some temper- 
ature Tbec < Tc into the unique single-particle state at 
the bottom of the spectrum. 



II. SPIN-ORBIT COUPLED BOSONS 

We consider a two- component Bose gas with Rashba 
spin-orbit coupling and a two-body interaction which is 
invariant under rotations around the z-axis in spin- space. 
The Hamiltonian is given hy H = Hq -\- Hint, with 
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where k±^ = k^x^kyif is the projection of the wave- vector 
k onto the x?/-plane, and the components of the vector 
operator cr are the usual Pauli matrices. The wave- vector 
/co measures the strength of the spin-orbit coupling. The 
invariance of the interaction under spin-rotations around 
the z-axis implies that the function /7cricr2(^i7 ^2 5 ^2, ^i) 
is symmetric with respect to the simultaneous permuta- 
tion of its incoming and outgoing labels^^ 

For simplicity, we assume that the bare interaction is 
momentum independent, so that the interaction is char- 
acterized by three different coupling constants 



^t = t/n(0); g^ = Uu{0); g^ = Un{0) = U^^{0), 



and Hint simplifies to 
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+9l(^l,i(^l^i<^k2iak,i + 2g_Lal^^al^^ak^^ak,t] ■ (2-5) 
Introducing the usual field operators 



Mr) = ^j:e^''- 



O^ko 



(2.6) 



and the corresponding density operators paix) 
tljl{r)ilja{r) the interaction can be written as follows, 

1 



Hi] 



d r g^: p^{r) : +^; : p^{r) : 



+2^± : p^{r)pi{r) : 



(2.7) 



where : . . . : denotes normal ordering. Assuming for sim- 
plicity g^ = gi = 9\\^ we may write the interaction as 



Hi] 



d^r 



9p ■■ P^{r) : +g^ : a'^{r) : 



(2.8) 



where p = p^ -\- pi and a = p^ — pi represent the density 
and the spin density, and the associated couplings are 



9p 



1 



1 



{g\\^g±)^ 9a = 7:{g\\-g±)- (2.9) 



From Eq. (2.8) it is clear that for ga < (i.e., g± > g\\) 
states with a finite spin-density are favored; an example is 
the standing wave spin-striped state proposed in Ref. [4J , 
where the bosons condense simultaneously into two mo- 
mentum states with opposite momenta on the low-energy 
manifold in momentum space. On the other hand, ga > 
favors states with vanishing spin density, such as plane 
wave condensate where the bosons condense into a single 
momentum state, or the charge striped states discussed in 
Refs. [6^ and 13j. In this work, we focus on the case where 
g± is negative. We show below that even for infinitesi- 
mally small g± < the system exhibits spontaneous fer- 
romagnetism in the plane of the spin-orbit coupling at 
sufficiently low temperatures. 

To set up our notation, let us review the diagonaliza- 
tion of Ho- Performing a momentum-dependent rotation 
in spin space around an axis 6k/\0k\ with angle \6k\, 
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we obtain 
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(2.12) 
We now choose the rotation matrix e^^^ such that it 
rotates the z-axis into the direction k± = k±/\k±\^ i.e.. 



k^ = e^^''z, 
which can be achieved by setting 



Ok 



-z X k^. 



(2.13) 
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Due to the rotational invariance of the scalar product, 
we may hence write 

k^ • {e^^'^cr) = {e^^^'z) • {e^^"" ct) = z ■ a = (j\ (2.15) 

Our rotation matrix in spin space is then explicitly given 
by 
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where in the last line we have set kx = cos (^k and ky = 
^vdipk' In the new basis (which we shall call the helicity 
basis) the non-interacting part of the Hamiltonian is 



Ho = 2_^ z_^ Ekxaj^^akx, 



(2.17) 
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FIG. 1. (Color online) Graph of the energy dispersions (2.18) 
of the spinor Bose gas with Rashba-type spin-orbit coupling. 
The minimum of the lower helicity branch is a circle in the 
plane /c^ = with radius /cq. The spacing between contours 
scales quartically. 



with energy dispersions 
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Here vq = /co/m and the helicity index A = ± labels the 
two branches of the dispersion. A graph of these dis- 
persions is shown in Fig. fl] Note that the energy Ei^_ 
assumes its minimum —kQ/{2m) on a circle of radius k^ 
in the x?/-plane, while £^fc,+ is non- negative and vanishes 
only at /c = 0. Due to the momentum-dependent rotation 
in spin-space, the interaction vertices in the helicity ba- 
sis acquire a momentum-dependence. For completeness 
we give the properly symmetrized expressions for these 
vertices in Appendix A. For our mean-field calculation it 
is more convenient to work in the original spin basis. 



III. MEAN-FIELD THEORY FOR 
TRANSVERSE FERROMAGETISM 

A. Derivation of the mean-field equations 

To study transverse ferromagnetism we add a uni- 
form magnetic field h^ in the xy-plane, so that the non- 
interacting part of our Hamiltonian is now given by 



^o = ^(4t'4i) 



k^ - 2kok_i • a 
2m. 



h± • a 






(3.1) 



For convenience we measure the magnetic field in units 
of energy. The system exhibits spontaneous ferromag- 
netism if the magnetization remains finite when h^ -^ 0. 
The spin-rotational invariance with respect to rotations 
around the z-axis is then spontaneously broken. While in 
the symmetric phase the self-energies T^aaf are diagonal 
in the spin-labels, in the symmetry broken phase there 
are finite off-diagonal components E^^ and S^^. Within 



the self-consistent Hartree-Fock approximation the self- 
energies are independent of momentum and frequency if 
we start from a momentum-independent bare interaction. 
The mean-field Hamiltonian is therefore of the form 



HMF=Ho + ^(a^^,a^^) 
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Within the self-consistent Hartree-Fock approximation 
the self-energies are 
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where we have introduced the densities 
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(3.4) 
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Here the expectation values should be evaluated with the 
grand canonical density matrix associated with the mean- 
field Hamiltonian (3.2). Keeping in mind that p^i = 
Mx + iMy can be expressed in terms of the components 
of the transverse magnetization Mj_ = MxX + Myij^ we 
see that the Cartesian components of the off-diagonal 
self-energies are proportional to the corresponding com- 
ponents of the magnetization, 



5]. = -(S^i + S;^)=^±M„ (3.6) 

It is convenient to define in addition the self-energies 

^z = 2 (^n - ^u) = 9tPt - 9iPi - ^ ' 



■-{Pt-Pi). (3.8) 



So = -(S^t + ^u) = 9tPt + 9iPi ^Y^^^^ ^^)' ^^'^^ 



and the wave-vector 



p= {h±- ^)/vo, 



(3.10) 



where T, — T^^^ ^^yit ^^z^ is proportional to the in- 
ternal magnetic field induced by the interaction. The 
Hamiltonian can now be diagonalized via a momentum- 
dependent rotation in spin space of the form (2.10). The 
rotation matrix can be written as 



g-^O-^fc ^ 
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FIG. 2. (Color online) Graph of the energy dispersions (3.13) 
of the spinor Bose gas with Rashba-type spin-orbit coupling 
and an effective magnetic field h± —Ti = vop pointing in the 
direction of the positive x-axis. The lower helicity branch has 
a unique minimum at ko = kop/\p\. The spacing between 
contours scales quartically. 



where the rotation angles in the presence of an effective 
magnetic field h± — X) = vop are now given by 
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The energy dispersions of the eigenmodes are 
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where A = ± labels again the helicity of the modes. These 
dispersions are shown graphically in Fig. |2j Obviously, 
for any finite p the degeneracy of the corresponding dis- 
persion without magnetic field shown in Fig. [T] is com- 
pletely removed, so that E^,- now has a unique minimum 
at ko = kop/\p\. 

To derive a self-consistency equation for the transverse 
magnetization, we simply evaluate the off-diagonal den- 
sity (3.5), using the grand canonical density matrix asso- 



ciated with the mean-field Hamiltonian on the right-hand 
side. We thus obtain 
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where 
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is the average occupation of the mode with energy E^x. 
Below we shall work at constant density, so that we 



should eliminate the chemical potential /i in favor of p. 



Therefore we need the diagonal densities (3.4), 
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implying that the total density is 
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To show that the phase transition to the ferromagnetic 
state is continuous, it is useful to calculate the grand 
canonical potential, which in a mean-field approximation 
is given by 

n{T,ii,h) = r^ ^ In [l - e-''(^--^)] - {-Hint), 

k \=± 

(3.18) 
where the expectation value of the interaction part of the 
Hamiltonian is 



{nint)/v 
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(3.19) 



and the second line holds for ^^ = ^^ = ^y . For simplic- 
ity, let us now choose h± = hx so that M± = Mx. To 
explore the possibility of spontaneous transverse magne- 
tization at constant density, we should consider the Gibbs 
potential 

G(T, p, M) = n{T, /i, h) + V{/j.p + hM), (3.20) 

where p = fi{p, M) and h = h{p^ M) should be deter- 
mined by inverting the equations 

' = -vJ-p^ ^ = -vml- 



(3.21) 



To study spontaneous transverse ferromagnetism, we 
shall later take the limit h ^ 0. For simplicity, we assume 
that g^ = gi = g\\, so that Pt = Pi = p/2 and So = 
(^'11 + g±/'^)p- As a consequence E^ = and Pz = 0. 
Our self-consistency equation (3.14) then reduces to the 



following equation for the transverse magnetization. 



2VJ-^ \k^^p\ 
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Note tha t p in the ri ght-han d side depends again on M± 
via Eqs. (3.10) and ( 3.6|3.7 ); for g^ = g^ and /i ^ the 
relation between p and M±_ is simply p = —g^M^/v^. 
To determine the order of the phase transition for /i ^ 0, 
it is sufficient to consider the change in the free energy 
F(T, p) = 1^ + Vpp for arbitrary magnetization M, 

AF(T, pM)= F{T, p,h = 0)m^o - F{T, p,h = 0)m=o. 

(3.23) 
The physical state of the system at vanishing external 
field is determined by dAF{T, p; M)/dM = 0, which 
i s ano ther way of deriving the self-consistency equation 
(3.22) for the order parameter M. 



B. Spectral densities 



To evaluate the integrals appearing in Eqs. (3.17) and 



(3.18), it is useful to introduce the density of states 



^A 



(^'^) = ^E^(^-^^a) = ;^E^(^-^^'-P'A), 



y^ . ... y ^^ 
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where p = {h — g±M)/vo and we have shifted k' = k -\-p 
on the right-hand side. The density equation ( 3.17| ) can 
then be written as 



/OO -j^ 



(3.25) 



while our expression (3.18) for the grand canonical po- 



tential per volume becomes 
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It is also useful to rewrite the self-consistency equation 
for M± = Mx in terms of a general ized s usceptibility as 
follows. Shifting k' = k -\- p in Eq. (3.22) we obtain 



M: 
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We now introduce cylindrical coordinates in A^'-space and 
perform a partial integration in the angular part. Rear- 
ranging terms we find that the self-consistency equation 
(3.27) can be written as 



h 



1 
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where the susceptibility X-l{M) is defined by 



(3.28) 



It is easy to see that x±(^) > 0? implying that, only 
for g± < 0, the magnetization M can remain finite for 
h ^ 0. From now on we shall therefore assume that g± 
is negative. In this case the magnetization in the broken 
symmetry phase satisfies the self-consistency equation 



X±{M) 



1 
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To evaluate x±{M), we introduce the weighted density 
of states, 

-^(^.^) = ^E^^(-^-.^a). (3.31) 



Then we may write 
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y_^de^AcT.(e,p)p(j3;^^-^.(3.32) 



Introducing cylindrical coordinates in /c'-space in the 
above integrals defining ux{e^p) and <jA(e,p), the inte- 
grations over k'^ and \k'j_\ can be carried out exactly so 
that we can write these functions as one-dimensional an- 
gular integrals. The results can be written in the scaling 
form 



z^A(e,p) . /e-So 
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In Appendix B we give explicit expressions for the di- 
mensionless scaling functions i>x{e^p) and ax{e^p) as one- 
dimensional integrals. In fact, for negative e the remain- 
ing angular integration in the expressions for z>_ (e, p) and 
( 7_ (e , p) can also be done analytically, see Eqs. ^B2[ ) and 
( |B4| ). Graphs of the scaling functions are shown in Fig.js] 

The behavior of the spectral densities V-{e^p) and 
(j-{e^p) associated with the negative energy branch is 
rather interesting. Both functions vanish if e is smaller 
than the lower threshold energy 



-eo + So - ^oP. 



(3.36) 



Recall that in the absence of an external magnetic field 
p = —T^x/vq = {—g±)M/vo. For energies slightly above 
the lower threshold we obtain from Eqs. (B5, B6), 
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where p = p/ko. Eqs. (3.37, 3.38) can also be derived by 
expanding the energy dispersion Ej^_ of the lower branch 
around the minimum ko = k^x to quadratic order. 
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This approximation is only accurate for |£^fcQ+q _ — e_ | <C 
2voP] the energy surface E^- = e can then be approxi- 
mated by an ellipsoid. However, for higher energies the 
topology of the energy surface changes, as illustrated in 
Fig. |4] With increasing energy the ellipsoid distorts into 
a bean-shaped surface until the two ends of the bean meet 
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FIG. 3. (Color online) (a) Graph of the scaling functions 
^±(^^V) of the density of states defined via Eq. (3.33) for p = 
0.1, see Eqs. (Bl |B2 ). (b) Graph of the scaling fun ctions 
5"± (e, p) of the weight ed densit y of states defined via Eq. ( 3.34 ) 
for p = 0.1, see Eqs. (B3 B4). The solid lines correspond to 

-1 branch while the dashed 
1 branch. The insets 



the scaling functions in the A = 

lines are the scaling functions in the A ■■ 

show a closeup of the negative energy part of the spectral 

functions in the lower helicity branch. We use e- = e-/eo 

and e* = e* /co- 



at a critical energy. For higher energies, a hole emerges 
in the energy surface so that it assumes the topology of 
a torus. At the critical energy e* where the topology of 
the energy surface changes the spectral densities have a 
cusp. From the exact expressions for the spectral densi- 
ties given in Appendix B it is easy to see that the critical 
energy is 



e* = 
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VoP = e- 



2voP- 



(3.40) 



A similar transition in the topology of the Fermi sur- 
face of metals as a function of external pressure has been 
discussed a long time ago by Lifshitz.^^ Close to such a 
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FIG. 4. (Color online) (a) Surface of constant energy ^fc,- = e 
of the lower helicity branch in the ferromagnetic phase just 
above the lower energy threshold e_ = — (1 +p)^. (b) Same 
as (a) but enlarged by a factor of seven, (c) Evolution of the 
constant energy surface for different energies. The scales are 
the same as in (a). All plots are for p = 0.1. 



transition, Lifshitz predicted anomalies in the thermody- 
namics and the kinetics of the electrons. Therefore we 
expect that in phases with spontaneous transverse ferro- 
magnetism the kinetics of spin-orbit coupled bosons with 
energies close to e* is rather unusual. 



C. Solution of the mean-field equations 



For the numerical solution of the above mean-field 
equations it is useful to introduce the dimensionless den- 
sity, magnetization, susceptibility, and interaction as fol- 



lows, 
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(3.41b) 

(3.41c) 
(3.41d) 

We also introduce the dimensionless energy u = {e — 
e_)/eo which is measured relative to the bottom of the 
lower helicity branch, and define 

(3.42) 

(Tx{uJ,p) = =CrA(-(l+p) +CJ,p). 



z^o 



(3.43) 



Finally, we introduce the dimensionless temperature 

T = T/eo, (3.44) 

and the fugacity 

^^g(/i-e-)/T^ (3.45) 



With this notation the density equation ( 3.17[ ) can be 
written as 



(3.46) 






while the self-consistency equation ( 3.3Q| ) for the dimen- 
sionless order parameter p = p/k^ = —g^M/2 becomes 



X±(r,2;,p) = -l/^±. 



with 



YrL ^^E^^^(^'^)l^;W^3^- (3.48) 



Finally, the dimensionless free energy /{r^p^p) = {Q 
liN)/{yvQe^) can be written as 



f{r,p,p)=T / da;^i/A(^,p)ln 
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(3.49) 



where p = /i/eo. Note that at constant density we should 
determine /i as a function of p and p. 

Let us first discuss the critical temperature Tc = Tc/cq 
below which the system exhibits spontaneous transverse 
ferromagnetism. According to Eq. (3.47), for a given den- 



sity p the critical temperature is determined by 

X±(^c,^c(rc,p),p = 0) = -l/g±, (3.50) 
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FIG. 5. (Color online) (a) Critical temperature for trans- 
verse ferromagnetism as a function of the dimensionless cou- 
pling constant \g\_\ for three different densities. The dots 
have been obtained from the numerical solution of the mean- 



(3.47) field equations (3.50) without further approximation, while 



the so lid lines represent the low-temperature approximation 
( |3.55 ). (b) Critical temperature as a function of density for 
different values of the interaction. 



where the fugacity Zcijc^p) at the critical point is deter- 
mined by 



coo ^ 



(3.51) 



Numerical results for Tc as a function of \g±\ for differ- 
ent densities are shown in Fig. Is] (a), while in Fig. Is] (b) 
we show the critical temperature as a function of den- 
sity for different values of |^^|. The numerical results in 
general are obtained without approximating the spectral 
densities and choosing g\\ such that Eq vanishes. Note 
that for small \g±\ the critical temperature approaches 
zero with infinite slope. In this regime it is easy to ob- 
tain an analytic expression for the critical temperature. 
Assuming that the temperature is small compared with 
eo (corresponding to r <C 1) we may neglect the contri- 



bution of the upper helicity branch and approximate the 
spectral densities by their leading asymptotics for fre- 
quencies close to the bottom of the lower helicity branch 
given in Eqs. (B13) and (B14). I n the symmetric phase 
where p = the density equation (3.46) then reduces to 



P 



duj — 

6 



h 



-rln(l-z). 



(3.52) 



We conclude that in the regime where the critical tem- 
perature Tc is small compared with unity, the critical fu- 
gacity is given by 



1 



-p/t 



(3.53) 



For Tc ^ p this is exponentially close to unity, which 
is a consequence of the finite density of states of the 
spin-orbit coupled Bose gas close to the bottom of the 
lower energy branch. To determine the critical tem- 
perature as a function of the densi ty, we calculate the 
transverse susceptibility from Eq. (3.48) with p = 0, 



using the approximation (B14) for the spectral density 
a_(a;, 0) = cr_(— 1 +0^,^ = 0) for frequencies close to the 
bottom of the lower helicity branch. 



X± 



4 70 



P^l^ 



duj- 



1 



o/r 



41 



oP/-^ 



(3.54) 



Hence, in the paramagnetic phase the transverse suscep- 
tibility becomes exponentially large at low temperatures. 
As a consequence, for any finite attractive interaction 
g± < 0, we ca n find a solution of the self-consistency 
equation ( 3.4^ ) at sufficiently low temperatures. Com- 
bining Eqs. (3.50) and (3.54) we obtain for the critical 
temperature 



P 



In (1 + 4/1^x1)' 



(3.55) 



From the derivation of this expression it is clear that 
Eq. (3.55) is valid as long as Tc <C 1, which can always 



be satisfied for sufficiently small densities. The approxi- 
mation (3.55) corresponds to the solid lines in Fig. [s] 

Next, let us discuss the low-temperature phase r < Tc 
with spontaneous transverse ferromagnetism. In Fig. |6] 
we show the change 



A/ = /(r,p,p)-/(r,p,0) 



(3.56) 



in the dimensionless free energy defined in Eq. (3.49) as 
a function of the order parameter M = 2p/\g±\ for three 
different temperatures. For r < Tc the free energy contin- 
uously develops two degenerate minima, corresponding 
to the Hartree-Fock solutions ±M. The phase transi- 
tion to the magnetic state is therefore second order. The 
transverse magnetization as a function of temperature for 
three different densities, obtained numerically, is shown 
in Fig. [7[ In the regime where Tc <C 1 the behavior of 
the magnetization for temperatures close to the critical 
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FIG. 6. (Color online) Dimensionless free energy A/, defined 
in equation ( |3.56[ ) , as a function of the dimensionless magneti- 
zation M for three different temperatures, r /tc = 0.9 (dotted 
line); t/tc = 0.95 (dashed line); t/tc = 1.1 (solid line). For 
all plots we have used |^±| = 5 and p —1. 
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FIG. 7. (Color online) Transverse magnetization as a function 
of temperature. The plots are for \g±_\ — 1. The dots are 
obtained numerically, while the solid lines correspond to the 
analytic result (3.57) which is only valid for Tc — r <^Tc ^1. 



temperature can be calculated analytically, as shown in 
Appendix C. We find 



M > 



ii^J)Mi + ^] 



Tr-T 



(3.57) 



The reason we obtain the usual mean-field exponent /3 = 
1/2 is of course related to the fact that our calculation is 
based on the Hartree-Fock approximation. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have shown that in the spinor Bose 
gas with Rashba-type spin-orbit coupling an arbitrarily 
weak attractive interaction between bosons with opposite 
spin triggers a ferromagnetic instability for temperatures 



below some finite temperature Tc if the density of the 
bosons is fixed. Note that spontaneous ferromagnetism 
in electronic systems usually appears only if the relevant 
interaction exceeds a finite threshold. ^^ The fact that in 
the spinor Bose gas such a threshold does not exist is 
related to the singularity of the Bose function for small 
energies in combination with the finite density of states 
in the non-magnetic phase due to the spin-orbit coupling 
of the Rashba-type. Because for T < Tc, the density of 
states exhibits the usual y^-behavior at low energies, at 
some temperature below Tc the bosons eventually con- 
dense into the single-particle state with the lowest mo- 
mentum ko, which is unique in the ferromagnetic phase. 
In the weak coupling regime, we may estimate the crit- 
ical temperature for BEC by using the critical tempera- 
ture for fr ee bosons with anisotropic dispersion given by 
Eq. KM, 



Tbec = 



27r 
m 

27r 
m 



C(3/2) 

P 

C(3/2) 



2/3 



2/3 



P 



ko 



1/3 



-Py 

\g±\M 



kl/m^\g^\M 



1/3 



(4.1) 



For consistency, we should require that Tbec < Tc, which 
is satisfied at weak coupling where p <C /cq. We thus 
conclude that the ground state of the spinor Bose gas 
with Rashba spin-orbit coupling and attractive interac- 
tion ^x < is a ferromagnetic Bose-Einstein condensate. 

In principle, an attractive interaction g_\_ < could 
also trigger an instability in the particle-particle channel, 
leading to a pair condensate at low temperatures. How- 
ever, as shown in Appendix D, at least for sufficiently low 
densities the ferromagnetic instability has a higher crit- 
ical temperature if the relevant coupling constants have 
the same order of magnitude. 

Let us also point out that the constant gj_ in our mean- 
field calculation should be considered as an effective 
low-energy interaction in the opposite-spin particle-hole 
channel. This interaction includes renormalization effects 
from high-energy fluctuations in all channels, so that it 
is in principle possible that the effective low-energy cou- 
pling is attractive even if we start from a repulsive bare 
interaction. Note that recently Gopalakrishnan et al.^ 
performed a momentum shell renormalization group cal- 
culation for the spin-orbit coupled Bose gas whose dis- 
persion exhibits a minimum on a circle in momentum 



space; at finite temperature they found evidence for an 
attractive renormalized coupling in the particle-particle 
channel, although the corresponding bare interaction is 
repulsive, suggesting an instability towards pair conden- 
sation. However, a proper renormalization group calcu- 
lation of the effective low-energy coupling of the spinor 
Bose gas with Rashba-type spin-orbit coupling, taking 
high-energy fluctuations in all channels and all relevant 
and marginal couplings consistently into account, still re- 
mains to be done. Given the fact that the system belongs 
to the Brazovskii universality class, the rescaling step in 
the renormalization group transformation is non-trivial 
and all two-body scattering processes describing particles 
with momenta on the low-energy manifold are marginal 
and should be retained. It should be advantageous to use 
the well deve loped m achinery of the functional renormal- 
ization groupP^I^^'^ to carry out such a calculation. 
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APPENDIX A: INTERACTIONS IN THE 
HELICITY BASIS 

In this appendix we explicitly give the interaction part 
of the boson Hamiltonian (2.2) in the helicity basis. For 



simplicity, we assume momentum independent bare in- 
teractions, see Eq. (2.4). To transform the interaction 



from the spin basis to the helicity basis, it is useful to 
introduce the notation a^ = ak- and bk = e~'^'^^afc+, so 



that our transformation (2.10) to the helicity basis can 
be written as 

akt = —^ [(^k- - e"^^'^afc+] = — = [a^ - hk] , (Al) 



V2 



v/2 



(^ki = —^ [e'^^^ak- + ak+\ = -^ 



[ak^hk]. (A2) 



Defining the coupling constants g^, g^ and g^ as in 
Eq. (2.4), we find that in the helicity basis the interaction 
part (2.5) of our Hamiltonian can be written as 
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T-Ci, 



V T. ^K 



V 



\ -\-k'^,k2+k 



h'-^h'^h2h\ 



, { i [rr°(l'. 2'; 2, l)a\,a\,a2a, + rf''(l', 2'; 2, l)h\,h\,h2hi 



+rf ''^(1', 2'; 2, l)at,at,;,26i + r^'^««(l', 2'; 2, 1)61,6^,0201 



+ 



2! 



rf °''(1', 2'; 2, l)a\,a\,a2hi + r^''''«(l', 2'; 2, l)6l,6^,62ai 
+r^«««(l', 2'; 2, l)6l,4,a2ai + Tf''\l\ 2'; 2, 1)4,6^,6261 
+rg'^^'^(l',2';2,l)al,6^,62ai} 
[7 X] '^fci+^,fe2+fci{7^n2^i('^i''^2;fe2,fci)[al'4'«2ai +6^,6^,6261] 



y 



k'.^k'^k2hi 



l(2!)2 



-— "2 /72(/ci, /C2; ^2, /ci)[a|,4,62^i + h{,h\,a2ai] 



^t 7.t 



Us{k[,k2] k2] ki)[a[,al,a2bi + 6{/62/62ai 



.t 7.t 



U4{k[;k2; /c2, /ci)[6{,a2/a2ai + ayh\,h2h 



-^Us{k[;k2;k2] ki)a\,bl,b2ai >. 
The properly symmetrized interaction vertices are 

rr^il', 2'; 2, 1) = rf''(l', 2'; 2, 1) ^ U,{k[,k',; k^, fci) 

2^2 ^ 4 ^ 



rg"''^(i', 2'; 2, 1) = rg^""(i', 2'; 2, 1 
rf "''(1', 2'; 2, 1) = C''"(i', 2'; 2, 1 
r^"""(i', 2'; 2, 1) = rf ^''(1', 2'; 2, 1 

rf'"'(l',2';2,l 



2^2 4 ^ 



/72(/ci,/C2;^2,^l) 

2^2 

U^{k[,k'2]k2]ki) 

_9l_ I ^ i((^^,+(^2/- 
2 ^ 2 

/74(/ci;/c2;^2,^i) 



+ e-i^2/)(e^^i +e^^2), 



■z^i, +e-*^2')(e*^i +e*^2)^ 



-'/?2-^i) ^ ^f^-i^i' _^ g-*'/'2' Ug*^i _ e^'^^A 



= U^{k[]k'2]k2]ki) 

_ 9^ j^ 9]^^-i{y^^,-\-ip2'-'P2-<Pl) _^ £i/g-*^l/ _ g-*'^2' Vg*^l _ g^^2\ 



r 



(A3) 

(A4a) 
(A4b) 
(A4c) 
(A4d) 
(A4e) 



in all angles. 

APPENDIX B: SPECTRAL DENSITIES 



In this appendix we discuss the dimensionless scaling 



The parametrization of the vertices in terms of five func- 
tions t/i , . . . , [/s is similar to the parametrization used 
by Ozawa and Baym;^ however, we find it convenient 
to introduce slightly different numerical prefactors in the 
second line of Eq. (A3) in order to simplify the combi- 
natorial factors due the permutation symmetries of the functions i>x{e,p) and ax{e^p) which are defined by writ- 
vertices in higher order calculations. Note that the in- ing the density of states and the w eighted density of 
teraction vertices are invariant under arbitrary rotations states in the scaling form (3.33, 3.34). The scaling func- 
around the z axis, corresponding to a shift cpk ^ ^k -^ o^ tion vxie^p) for the density of states can be written as 

I 



vx{e,p) = 9(6) I ^ + r J(pcos(^ - A) 



1 1 f pcosip — A 
— \ — arctan — 

2 TT \ ^l 



r27r 



27r^ 



+ 6(— e) / ^^(pcos(/? — A)6(pcos(/:) — A - 
Jo 



^~). 



(Bl) 
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For negative e the angular integration can be done analytically. For the lower energy branch (A = — 1) we obtain 



z>_(e,p) = 6(1 - V^ - p) + e{p - |1 - V^l) 



1 



e-1 



+ V^' - (1 



2\2 



e) 



, for e < 0. (B2) 



The scaling function for the weighted density of states d\{e^p)^ can be written as 



a,{e,p) = e(?) r S^ I [3(pcos^ - A)^ + e] 



e(-e) / 



27r 2 
^^ d(^ sin^ (/:) 



1 1 f pcoscp — A 
— \ — arctan — 

2 TT y y'g 



(pcos(/p — A) 

TT 



27r 2 



[3{pcos(p — A)^ + e] Q{pcos(p — A — a/^). 



(B3) 



For negative e the angular integration can again be performed analytically. For the lower energy branch (A = —1) we 
obtain 



(B4) 






3 + e + -p arccos 



, for e < 0. 



r 



Plots of the above scaling functions are shown in Fig. |3] 
The asymptotic behavior slightly above the lower thresh- 
old e_ = -(1+p)^ is, 



z>_(e,p) 



cr-{e,p) 



\/l+p 


/e+(l+p)2 




/ P 
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(B5) 



(B6) 



To discuss the thermodynamics close to the critical point, 
we need the expansion of the density of states and the 
weighted density of states for small p = p/ko to order p'^. 
We obtain 



^a(€,p) _ ~(o) / e 
where 



eo 



^0 \ , ~2~(2) / e 



eo 



(B7) 
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^ ^ -^ 2 7r(l + e 



2^2 



e(-e) 



Sx,-i5'{l - V^e), 



(B8) 



(B9) 



and S'{x) = j^S{x) is the derivative of the (5-function 
with respect to its argument. The order parameter ex- 
pansion of the weighted density of states is 



^a(€,p) ^ ~(o) /e 



^0 \ , ~2~(2) / e 
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(BIO) 
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(B12) 



At low temperatures only the leading asymptotics of 
these functions close to the bottom of the lower energy 
branch is relevant. Shifting the dimensionless energy as 
e = — 1 + w, we may approximate for |u;| ^ 1, 



i>^°\-i + uj)^e{uj), 



<tL°\-1+^) 



e{uj) 



(B13) 
(B14) 



i>^y{-l+io)!^h'{l-VT^^)!^5'{uj), (B15) 



<ji'\-l+w) 



1 

16 



[-e{uj) + 25{uj)+45'{io)], (B16) 



where we have used (5' (aw) = a ^S'{ui) and uj5'{ui) = 



APPENDIX C: ORDER PARAMETER 
EXPANSION 

For temperatures slightly below the critical tempera- 
ture the order parameter p = \g±\M/2 is small so that 



12 



we may expand all quantities to second order in pow- 
ers of p. Having expressed all quantities in terms of the 
spectral densities v\{e^p) and <jA(e,p), we use the order 
parameter expansion of these, given in Eqs. ( B7[ BIO), to 
generate the corresponding expansion of any other quan- 
tity. Using the density equation (3.46) we obtain for the 
chemical potential 



/i = /iO+M2P^ +C>(P'^), 



(CI) 



where the chemical potential for vanishing order param- 
eter is determined by 



^^X]4°^(^)^o(^)' 



with 



^o(e) = 



3(e+So-Ao)/T _ \ 



(C2) 



(C3) 



Here Hq = ^o/eo and JIq = /io/eo- The leading correction 
to the chemical potential in the symmetry broken phase 
is 



M2 M2 



/-oo^^Ea^a (^)^o(e) 



(2). 



^ I-oo d^Ex ^r(^~)^o(e)[no(6) + 1] 



(C4) 



The corresponding expansion of the susceptibility x± de- 
fined in Eq. (|3.48|) is 



~(0) , ~2~(2) 



o(/), 



(C5) 



with 



1 r'^ 

xT = -^J ^ de^ Aaf (e)no(e1[no(e~) + 1], (C6) 



and 



1 f^^ ( 

xf =-^y ^de^AJaf (e)no(?)[no(e~) + l] 

A 

_ ^^af\e)noCe)[noie) + l][2no(e~) + 1]|. (C7) 

And finally, the expansion of the dimensionless free en- 
ergy defined in Eq. (3.49) is 



n^i^N 



with 
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/CX) 
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fJ-oP, 
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T d~eY,vf\e)\u 
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-(e+i:o-Ao)/T 

(CIO) 



In the low-temperature regime r <C 1 the above ex- 
pressions can be evaluated analytically because it is then 
allowed to substitute the leading asymptotics of the spec- 
tral densities for ener gies close t o the bottom of the lower 
branch given in Eqs. (B13-B16). The zeroth order chem- 



ical potential and the corresponding fugacity zq are then 
given by 



zo 



,(Ao-So + 1)/t 



1 



-p/t 



while the correction to the chemical potential is 

1 e^/^ 



M2 



r(l - ^o) 



(Cll) 



(C12) 



r.(0) 



The zeroth order susceptibility x\_ is given by the right- 



hand side of Eq. (3.54), while the leading correction is 

2;o(3 + zo) 



~(2) 

xy 



8r2(l-zo; 



3 • 



(C13) 



The order parameter close to the critical point shows the 
usual mean- field behavior, 



p- (l-^o)i 



2pTc 



Tr-r 



zo{3^zo) 
At the critical temperature we may write 



I- zo 



-p/rc 



1 



4 



(C14) 



(C15) 



so that we arrive at Eq. (3.57) for the transverse magne- 



tization M = 2p/\g_\_\. Finally, the leading coefficient in 
the expansion of the free energy is 



/o = -r'^U2{zo) 



1 



^11 



91. 
2 



P +MoP, 



(C16) 



where Li2(^o) is the poly logarithm. The coefficient of the 
term proportional to p'^ is 



/2 



4 

9± 



zo 



1 



Zo 



1 

1^ 



~(0) 

XI 



(C17) 



Note that for r < Tc the right-hand side of this expression 
is negative, to that in the ferromagnetic phase the system 
indeed gains energy. 



APPENDIX D: PAIR CONDENSATION 

For attractive g± the system can also exhibit a pairing 
instability which competes with the ferromagnetic insta- 
bility discussed in the main text. In this appendix we 
show however that, at least for sufficiently low densities, 
the ferromagnetic instability is dominant. 
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If we decouple the interaction in the particle-particle 
channel, we obtain the mean-field Hamiltonian 

"Hmf - /uA/" = ^ [{Ek- - l^)alak + {Ek+ - n)blbk 



+Xoe''^''a-kbk + Xoe'^'^'^bla-k 



V 



IXol^ (Dl) 



In general, the boson-pairing order parameter xo = 
Ixok*" is complex, but all phases can be eliminated by 



setting ttfc 



^ah- The Hamiltonian can then be 



diagonalized by means of a Bogoliubov transformation, 

(D2) 
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The corresponding grand canonical potential is 
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The self-consistency equation for the order parameter is 



1 _ 1 ^ 1 



g/3(u;fc+A'Uo|fci|) — 1 2 

while the density equation is 
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(D7) 

The integrals are ultraviolet divergent and must be reg- 
ularized. One possibility is to eliminate the bare interac- 
tion in favor of the two-body t-matrix t± , which can be 
defined via 



1 1 



t± 9i_ V ^ 2efc 






(D8) 



The regularized gap equation is then 

^^^^\ek ojk 
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i\r /-^ i,u. 
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At the critical point we set xo = 0- At low temperatures, 
we may replace ^fc ~ cq — /i = 2eo in the prefactor on the 
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FIG. 8. (Color online) Critical density pc for pair conden- 
sation as a function of the two-body t-matrix |t±| (note that 
t± < 0, so —t± = \t±\) for fixed temperatures r. The sol id lin e 
depicts pc for vanishing temperature as given in Eq. (DIG). 



The dots repre sent numerical solutions of the regularized gap 
equation (D9) for the indicated temperatures. The shaded 



area underneath the solid line represents the regime where 
pairing cannot occur at any temperature. 



right-hand side. Then we obtain for the critical density 
Pc = Pc/i^o^o) at zero temperature 



Pc = 4 



1 



-{I'otA 



(DIO) 



For densities smaller than pc there is no pairing insta- 
bility at zero temperature. For finite temperatures we 
have solved the regularized gap equation ( |D9| ) numeri- 
cally, see Fig. [SJ Because the zero-temperature result 
(DIO) for the critical density forms a lower bound for the 



critical density at finite temperatures, we conclude that 
for sufficiently small densities, pairing cannot occur at 
any temperature. For densities smaller than this thresh- 
old the ferromagnetic instability discussed in the main 
text is dominant. 
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